Effect of adhesion geometry and rigidity on cellular force distributions 
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The behaviour and fate of tissue cells is controlled by the rigidity and geometry of their adhesive 
environment, possibly through forces localized to sites of adhesion. We introduce a mechanical 
model that predicts cellular force distributions for cells adhering to adhesive patterns with different 
geometries and rigidities. For continuous adhesion along a closed contour, forces are predicted to be 
localized to the corners. For discrete sites of adhesion, the model predicts the forces to be mainly 
determined by the lateral pull of the cell contour. With increasing distance between two neighboring 
sites of adhesion, the adhesion force increases because cell shape results in steeper pulling directions. 
Softer substrates result in smaller forces. Our predictions agree well with experimental force patterns 
measured on pillar assays. 

PACS numbers: 87.10.+e,87.17.Rt,68.03.Cd 



Adherent tissue cells react very sensitively to the physi- 
cal properties of their environment, including mechanical 
stiffness and the spatial distribution of adhesive cues [l[ . 
On fiat substrates, mechanical stiffness and adhesive ge- 
ometry can be controlled by using soft elastic substrates 
and microcontact printing of adhesive islands Q, re- 
spectively. Rigidity and geometry can be altered simulta- 
neously combining the above techniques [J, [5| or by using 
bio-functionalized pillar assays, in which cells adhere to 
the tops of an array of flexible micro-needles 0,0. Such 
biophysical approaches have revealed that stiffness and 
geometry sensing are closely related as they both involve 
forces being localized at discrete adhesion sites. 

Despite the tremendous experimental progress in this 
field, our theoretical understanding of the relation be- 
tween stiffness sensing, geometry sensing and cellular 
force distributions is still very limited. Here we analyze a 
simple theoretical model that describes the relation be- 
tween cell shape and force distribution. The shape of 
tissue cells adhering to discrete sites of adhesion is dom- 
inated by the formation of inward curved circular arcs of 
the non- adherent contour 

@, 0, EH E2- This shape fea- 
ture follows from a modified Laplace law of competing 
effective surface tension and effective line tension. Based 
on these concepts we calculate traction forces of station- 
ary cells as a function of adhesive geometry and stiffness. 




Model. Fully spread tissue cells typically flatten with 
only the nucleus sticking out, as shown schematically in 
Fig. Qji. Often they adhere to the substrate at sites of ad- 
hesion distributed along the cell boundary. In this case it 
is appropriate to consider an effectively two-dimensional 
(2D) model which parameterizes cell shape by its 2D con- 
tour f(l) with the internal coordinate /. After cell spread- 
ing is completed the cellular forces are mainly contrac- 
tile. Fig. [TJb depicts the physical forces acting at the 
cell boundary [8(. First the cell contour is drawn in- 



FIG. 1: (a) Cartoon of a typical tissue cell adhering to a flat 
substrate at four discrete sites of adhesion. In (b) we focus 
on the framed region, (b) Physical forces acting at the free 
boundary between two neighboring sites of adhesion: while 
the surface tension a pulls the contour inward, the line tension 
A pulls the contour straight. The black line represents an actin 
cable reinforcing the contour, (c) Here we consider three cases 
of increasing complexity: cells adhering to one large adhesive 
island, cells adhering to discrete sites of adhesion on a rigid 
substrate, and cells adhering to discrete sites of adhesion on 
a soft substrate. 



ward mainly due to spatially distributed tension in the 
actin cytoskeleton and the plasma membrane. In the 2D 
model, this effect is described by an effective surface ten- 
sion er, because the main effect is to reduce 2D surface 
area. Second the cell contour resists the inward pull. 
Often the cell contour is reinforced by the assembly of 
contractile actin cables connecting neighboring sites of 
adhesion. In the 2D model, this effect is described by 
an effective line tension A which contracts the contour 
between adhesion sites. In mechanical equilibrium the 
contractile forces exerted by cells are primarily balanced 
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follows from Eq. Q] 



FIG. 2: (a) Dimensionless force F/2X at equally spaced 
adhesion points as a function of the dimensionless span- 
ning distance /3 = ad/2\ for different opening angles <j> = 
7r/4, 7r/3, 7r/2, 27r/3, 37r/4, 7r from top to bottom, (b) Dimen- 
sionless force F/2X as a function of the dimensionless in- 
verse stiffness 7 = 2a /K for adhesion to regular polygons 
with <j} — 7r/4, tt/3, 7r/2, 2tt/3, 37r/4 from top to bottom and 
/3 = 0.1. 



on the adhesive substrate [12J . To derive the force distri- 
bution for a given adhesive geometry we therefore have 
to minimize the effective energy functional for the cell 
contour under the appropriate adhesion constraint: 



E= adA+ Xdl+ f ad ■ (r- r )dl 



(1) 



The first integral extends over the cell surface area A 
and the second and third along the cell contour I. The 
third term introduces the adhesion tension f a d as a La- 
grangian parameter for the constraint set by the adhesion 
geometry and described by r*o. Based on this model we 
calculate the cellular adhesion forces by solving the cor- 
responding Euler-Lagrange equation for the three cases 
depicted in Fig. [Tb. 

Continuous adhesion. For cells that adhere continu- 
ously to the substrate, the cellular shape f is fixed to 
r*o. The adhesion tension acting along the contour then 



fad = -[a+ \n]n 



(2) 



where n is the normal vector of the contour curve and 
k its local curvature. For convex shapes (k > 0), sur- 
face and line tension conspire to pull the cell contour 
inward, while for concave shapes (k < 0), the line ten- 
sion opposes the surface tension, thereby decreasing the 
adhesion force. For example, a cell adhering to a circu- 
lar patch of radius R applies a constant and inwardly di- 
rected adhesive tension a+X/R along the contour. Along 
straight boundaries (k = 0), e.g. along the straight lines 
of polygonal cells, the contribution of the line tension A 
vanishes and the adhesion tension is simply a. The con- 
tribution of the line tension localizes to the corners of the 
polygon. Approximating a corner with opening angle <j) 
by an arc with radius e and then taking the limit to a 
sharp corner by e — > 0, we can calculate the adhesion 
force F acting in the corner: 



F = lim 



(a + -)n(0)ed0 = 2Acos (|)n 6 (3) 
: e V 2 / 



where if = tt — <fi and n& points in the direction of the bi- 
secting line. The model predicts that the smaller </>, the 
larger the force pulling on the corner, with a maximal 
value of 2A when both arcs pull in the same direction. 
When considering a finite radius r of the adhesion, one 
has to set e = r rather than taking the limit e — > 0, result- 
ing in an additional contribution by the surface tension 
(7. Experimentally cells have been forced into various 
shapes and forces have been measured by using adhesive 
islands on compliant substrates 0] ■ On polygonal islands 
strong traction forces are measured at the corners. Our 
model qualitatively predicts this corner effect and makes 
quantitative predictions on the scaling with <p which can 
be experimentally tested in the future. 

Discrete adhesion. For cells that adhere to the sub- 
strate at discrete points of adhesion, the cell contour be- 
tween adhesion points is free and the shape equation re- 
sulting from Eq. [T] predicts the formation of circular arcs 
with curvature k = o~/A. The adhesion force follows from 
Eq. [3] by replacing the opening angle (f> with the actual 
pulling angle (f>* 7^ <p spanned by the two contour arcs 
pulling on the adhesion site. <fr* has to be derived from 
the shape equation. For three equally spaced adhesion 
sites with distance d spanning an angle <f> we can derive 
an explicit equation for the resultant force on the central 
adhesion site as a function of the adhesion geometry: 



F = 2A 



/JsinQ+v/i— ^ 



cos 



rib 



(4) 



with f3 — ad/2X being a dimensionless measure for the 
strength of the inward pull. can also be interpreted 
as dimensionless spanning distance d. Hence, the force 



scales again with the line tension A but now also depends 
on the spanning distance d and the surface tension a. 
Fig. [2^ demonstrates that the larger the spanning dis- 
tance d and the more acute the opening angle <f>, the 
steeper the inward pull and the closer the force comes to 
its maximal value 2A. This maximal value is reached at a 
critical pulling strength (3 C = sin(0/2) when the pulling 
direction has become so steep that the two arcs physi- 
cally touch each other (</>* = 0), at which point other 
physical processes will take over (e.g. pearling of tubular 
extensions @). For equally spaced adhesion sites along 
a straight line cf) = tt the adhesion force scales exactly 
linear with spanning distance d. As the opening angle 
(j) increases beyond 7r, the traction force may become a 
pushing force. For unequal spacing d\ 7^ e?2 between 
adhesion sites, a generalization of Eq. 0] can be derived, 
which shows that the force direction is now tilted towards 
the side with larger spanning distance. 

Elastic substrate. To study how in our model elastic- 
ity affects the adhesion force, we now consider discrete 
adhesion sites that move in a harmonic potential with 
spring constant K. For simplicity, we only consider reg- 
ular polygons, for which we find: 



F = 2A 



/3 70 + cos (I) y 7 2 +C0S 2 
7 ^ + cos2(|) 



n b (5) 



with 70 = sin(0/2) + 7Cos(</>/2). Here 7 = 2a /K is a di- 
mensionless inverse stiffness. Fig. [2)3 shows that the max- 
imal adhesion force is reached on a rigid substrate, 7 = 0, 
where Eq. [5] reduces to Eq. 2J As the substrate becomes 
softer, F decreases because the moving adhesion points 
effectively decrease the spanning distance d. Indeed it is 
well known that cellular traction is weaker on softer sub- 
strates 0. Like for rigid substrates the force decreases 
with increasing opening angle. Both these trends are 
also evident in Fig. [3h.-c, where we show computed cell 
shapes for triangles, squares and pentagons on stiff (red) 
and elastic (blue) substrates. 

Tension- elasticity model. Up to now, our results were 
derived with the assumption of constant line tension A. 
Recently it was suggested that this quantity has an elastic 
origin, i.e. A = EAu where u = (L — Lq)/Lq is the strain 
induced in the elastic contour with rigidity FA. Lo = ad 
is the resting length assumed to be proportional to the 
spanning distance d with a dimensionless resting length 
parameter a. The arc contour length follows from geo- 
metrical considerations as L = 2R a,rcs'm(d/2R), where R 
is the arc radius. Hence, the line tension A itself becomes 
a function of the adhesion geometry: 

d 



X(d) = EA 



2R 

— - arcsm , 
ad V 2R 



1 



(6) 




which increases with d. Thus also the radius of curvature 
R = X(d)/a increases with d, as indeed observed experi- 



FIG. 3: (a-c) Predicted cell shapes and adhesion forces for 
various adhesion geometries on stiff 7 = (red) and compliant 
substrates 7 = 1/15 (blue). Parameters /3 = 1/6 and 4> = 
7r/3,7r/2 and 3/5n from a to c. (d) Predicted cell shape and 
adhesion forces on a rigid substrate using constant line tension 
(red) and the tension elasticity model (green, a = 1, If = 20a) 
for an adhesion geometry with different spanning distances, 
namely d = a at the sides and d = y/2a at the top. In the 
constant tension model A was adjusted to result in the same 
force for d — a. 



mentally [ll|] . The distance dependence of the adhesion 
force now becomes a non-trivial function of the adhesion 
distance d. On the one hand, the line tension increases 
with d, thereby increasing the individual forces pulling 
on the adhesions. On the other hand, the arc curva- 
ture k = a/X(d) decreases with increasing d, leading to 
less steep pulling directions and therefore reduced overall 
force. In Fig. [3]i we compare the results of the tension 
elasticity model (green) to the constant tension model 
(red), where all arcs have the same curvature k = a/X. 
For the tension-elasticity model, the two top arcs across 
the diagonal have a reduced curvature K, because here the 
spanning distance d is increased by a factor of \ 2. The 
total force magnitude slightly decreases due to less steep 
pulling directions, although the tension in the individual 
arcs actually increases. 

Extracting model parameters from pillar assays. In 
general one expects the surface tension a to be a global 
quantity and the line tensions A to be different in different 
arcs. Our model suggests a simple procedure to estimate 
numerical values from the shape geometry. Fig.[4k, shows 
experimental data for an endothelial cell on a pillar array 
(Fig. 6 fromRef. 0). The contour was fitted by circular 
arcs using the procedure from Ref. Because here we 
focus on contour effects, we exclude all arcs which might 
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FIG. 4: (a) Predicted contour forces (red) and measured pil- 
lar forces (white) for an endothelial cell cultured on a pil- 
lar array [g]. Force direction is determined by cell shape 
geometry and force magnitude scales with a as determined 
from a least-square fit. Our procedure allows to predict the 
forces within each arc A = (8, 9, 15, 25, 16, 9, 12, 13, 38, 19) nN 
along the contour (counter-clockwise starting at the asterisk). 
Scales: Pillar spacing 9 /jm; force with plus 14 nN. (b) Force 
directions of pillar forces (black) and predicted contour forces 
(red). 



interact with other arcs due to close proximity or which 
might be distorted by internal stress fibers. The force on 
a pillar is the vector sum of the two adjacent arc forces. 
Assuming constant surface tension, the ratio of arc radii 
equals the ratio of line tensions, R1/R2 = Ai/A 2 , and the 
respective pulling angles /3j are given by cos (3i — Ri/2di. 
Thus, the resultant directions of the contour forces F are 
determined by R and d only. Fig. [4)d shows that we ob- 
tain excellent agreement between the predicted directions 
of contour forces and the measured directions of the pil- 
lar forces (error bars result from an estimated 10 percent 
uncertainty) . Because the surface tension determines the 
magnitudes of the contour forces, we extract a value of 



a ~ 2 nN///rri by a least-square fit to the experimental 
data, which is at the upper limit of the range of values re- 
ported earlier for cortical tension. The predicted contour 
forces (red) are shown in Fig. [4^ and compare favourably 
with the measured pillar forces (white). 

Discussion. Our model shows that shape and forces 
are closely related in cell adhesion and demonstrates how 
the spatial distribution of adhesion sites determines the 
forces acting at sites of adhesion. The effect of these 
forces might be further amplified by feeding into force or 
displacement-dependent regulation of cytoskeleton and 
adhesion sites, which should be included in future mod- 
elling. For example, Eq. [5] predicts that displacement 
F/K first increases linearly and then saturates as a func- 
tion of inverse stiffness, with important consequences for 
strain homeostasis. In this way, our model can be used 
to derive quantitative predictions how cell behaviour and 
fate can be steered by adhesion geometry and stiffness. 
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